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TECHNICAL NOTE 3289 

THE MINIMIZATION OF WAVE DRAG FOR WINGS AND 
BODIES WITH GIVEN BASE AREA OR VOLUME 
By Max. A. Heaslet 


SUMMARY 


The minimization of wave drag for thin aerodyn ami c shapes carrying 
no lift is studied for conditions under which either hase area or volume 
is specified. When volume alone is given, the analysis is limited to 
shapes with straight trailing edges normal to the stream direction. The 
problem of m inim ization is reduced to one of finding a two-dimensional 
harmonic function with known "boundary conditions. 

In several examples the theory is applied to the calculation of 
minimum drag. For given hase area, general formulas are found that cover 
as special cases quasi- cylindrical "bodies of revolution, w ing s having 
plan forms with fore and aft symmetry, slender Bodies, and certain classes 
of yawed wings. The drag can in fact he determined from a unidimensional 
flow analysis in a duct of known shape. For given volume, minimization 
of the external wave drag of a ducted "body of revolution of arbitrary 
radius is achieved in closed analytic form. 

In two cases, the determination of surface shape corresponding to 
the minimum drag is carried to completion. 


INTRODUCTION 


To seek conditions under which the wave drag of a given wing or "body 
is m inimi zed is to seek conditions for economical supersonic flight. It 
is also a common experience, in the study of such problems, to find t hat 
a gratuitous economy appears to affect the analysis itself. Almost invar- 
iably, simplicity characterizes the final forms of the results in compari- 
son with predictions carried out for wings and bodies chosen with less 
discrimination. In the present paper, the minimization of wave drag for 
aerodynamic shapes carrying no lift is studied. Conditions must, of 
course, be fixed in order to prescribe the problem and, from a practical 
point of view, this choice is less obvious than in studies of lifting 
configurations where a given weight is to be supported aerodynamically. 

The cases treated here apply the conventional constraints on base area 
and enclosed volume to a variety of shapes. For a large class of wings 
and bodies, the above-mentioned simplicity is especially apparent in the 
case of given base area for , as will be shown, the general expression 
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for minimi mi drag assumes tlie most elementary form possible while at the 
same time ret aining the relevant parameters and being dimensionally 
correct. 

The starting point of the present work is the expression for drag 
given by G. N. Ward (ref. 1) in his study of thin lifting bodies, that 
is, wings and bodies for which lin earized supersonic flow theory applies. 
The body shape is assumed to be enclosed by a characteristic surface gen- 
erated as the envelope of both the downstream- facing Mach cones, with 
vertices on the forward edge of the body, and the upstream-facing Mach 
cones, with vertices on the trailing edge of the body. Wave drag (plus 
vortex drag when lift is present) is then given by a control- surface 
integral of the induced velocities over the downstream portion of the 
Mach envelope. This particular control surface has analytical advantages 
similar to those exploited by R. T. Jones (refs. 2 and 3) in the use of 
combined flow fields. Jones adopts a perturbation potential equal to 
the sum of the potentials in forward and reverse flow. He then shows, 
for example, that the necessary condition for minimum wave drag is, for 
a plan form of given base area, that the pressure in the combined flow 
field be a constant over the plan form. It follows that locally the 
combined-flow potential is a two-dimensional harmonic function. If, how- 
ever, the entire analysis is committed to the use of combined flows, 
details of body shape (or, in the lifting case, surface loading) are lost. 
Along the Mach envelope used by Ward the perturbation potential in forward 
flow is equal in magnitude to its value in the combined flow field. Drag 
minimization then determines conditions on the control surface and for the 
problem of given base area the potential on the surface differs from a 
harmonic function by a known amount. The conventional perturbation poten- 
tial is retained but the determination of the body shape is still not 
direct. Mathematically, one needs to invert an integral equation and the 
question as to uniqueness of solution arises. The final examples in the 
present paper will be concerned with the construction of minimum-drag 
bodies of revolution from the known conditions on the control surface. 

In these special cases the inversions of the integral equations are easily 
carried out. 


ANALYSIS 


The following, analysis is divided into two sections. First, integral 
relations are derived that determ in e the drag, base area, and volume of a 
thin body in terms of induced velocities on the rear Mach envelope. Sec- 
ond, variational methods are applied so as to minimize wave drag subject 
to the imposed constraints. The perturbation potential on the Mach enve- 
lope is in this way related to a two-dimensional harmonic function 
satisfying known, boundary conditions. 
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Relations for Drag, Base Area, and Volume 


It is assumed that the given body 
deviates slightly from a planar or a 
cylindrical reference surface with ele- 
ments passing through the leading edge 
of the body and extending back parallel 
to the stream direction. As shown in 
sketch (a) , the body reference surface 
is denoted £ 0 , the over-all length 
is l, and a Cartesian coordinate 
system is to be used with the origin 
fixed at the foremost point of the 
body and the x axis alined with the 
free-stream direction. The constant 
free- stream velocity, Mach number, and 
density are U^,, 1^, and p^, respec- 
tively. The perturbation velocity components u, v, w induced by the 
body in the x, y, z directions are given by the gradients of the per- 
turbation velocity potential <p(x,y,z), that is, by cp x (x,y,z), qr(x,y,z), 
qp z (x,y,z) . Since supersonic small-disturbance theory is assumed to apply, 
the flow field is governed by the linearized differential equation 

P < Pxx “ c Pyy ” ‘Pzz = 0 

where 0 2 = Mm 2 - 1. 

Three types of reference surfaces occur most often in practice: 
first, planes which are associated with the study of wings; second, cir- 
cular cylinders which lead to the study of internal or external flow 
around quasi-cylindrical bodies; third, lines parallel to the stream 
direction. The last case is associated with slender-body theory and is 
merely a limiting form of the previous cases. The thickness distribution 
of the wing or body is fixed by boundary conditions prescribed on the ref- 
erence surface. In order to avoid difficulties conce rnin g gaps or holes 
in the body surface it will be assumed that unique leading and trailing 
edges exist and that the thickness distribution does not vanish between 
these extremities. 

In sketch (a) the characteristic surfaces enclosing the body are 
indicated. The front portion, E l7 is the envelope of the Mach cones 
stemming back from the leading edge of the body and the rear portion, 
is the envelope of the Mach cones facing forward from the trailing edge. 
The surfaces Ii and S 2 intersect along the space curve Ii. As shown 
by Ward (ref. 1) , the wave drag can be expressed in terms of an integral 
over the surface E 2 . This follows from an application Of momentum prin- 
ciples to the three-d im ensional region bounded by E 0 , Ex, and E 2 . The 
force on the body is expressed vectoria l! y in the form (see, e.g., ref. 4, 

p. 222) 



So 


Sketch (a) 
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■where the subscript oo denotes free-stream conditions, p and p are local 
static pressure and density, v is the local perturbation velocity vector, 
and the surface integrations extend over the "bounding surfaces. In small- 
disturbance theory the approximate relations 


p/p ro = 1 - *C(u/Uj 


and 


P-P Q 


“Poo 


U m u + i (_p 2 u 2 +v 2 +-vr 2 ) 


may be used. 


If cos(v,x), cos(v,y), cos(v,z) are the direction cosines of the 
inner normal v to the enclosing surface, drag (wave plus vortex drag) 
is, to the order of the approximations. 



D = ~ -vr // 1 (vS4-v2) + u 


0 2 u-!* 


_2w cos ( v ; z ) 

cos(v,x) cos(v,x) 


]} 


cos(v,x) dE 


( 2 ) 


An essential simplification of the drag formula follows if one 
introduces the function X(y, z) where 

X. = cp[f(y,z),y,z] 

and x = f(y,z) is the equation of the characteristic surface Z 2 . The 
function X is thus the value of the perturbation potential on Z 2 . The 
direction numbers of the normal on this surface are given by 

cos(v,x) :cos(v,y) :cos(v,z) = -l:f y :f z (3) 


and the relation 


holds . Since 


fy 2 + f Z 2 - P 2 


Xy = V + ufy x z = w + uf z 


the drag formula becomes 
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JJ (>y 2 + X z 2 )cos(v,x)dS 




(ha.) 

(to) 


The latter relation is a direct consequence 
of Green’s theorem in the two-dimensional 
cross-plane; V 2 is the Laplacian operator 
(d 2 /d y^)+(8 2 /8z 2 ) and n is the inner nor- 
mal to the boundary curves in the yz 
plane. As indicated in sketch (b) , the 
integration in equation (4b) extends over 
the area S which is the projection of S 2 
on the plane x = const. Along the outer 
curve C x , which is the trace of the curve 
T x in the same plane , the function X 
vanishes since also lies on the Mach 

cones from the leading edge. The inner 
curve C 2 is the trace of the cylindrical 
reference surface. 



y 


In order to evaluate the base area external to the reference sur- 
face £ 0 , it is sufficient to apply over the surfaces Z Q) Si, and Z z 
the first-order form of the integral expression for continuity of mass 
flow. This conservation relation is, for compressible flow. 



(U«* + cp)<a£ = 0 

dv 


and to the order of the present linearized analysis becomes 



v cos(v,y) w eos(v,z) 

u co cos (v ,x) ho, cos(v,x) 


cos(v,x)dZ = 0 


where the two integrations extend, as indicated, over So and. the rear 
Mach envelope Z 2 , respectively. In linear theory the term (l/Uco) (dcp/dv) , 
appearing as the first integrand, is equal to the slope, relative to the 
stream direction, of the body surface. The x-wise integration of the 
first integral thus yields the difference in body ordinates at the trail- 
ing and leading edges and the complete integration can then be written 
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2o 

where A is the incr ement of area between the nose and tail of the "body. 
In the notation of equations (3) and (4) , continuity of mass flow thus 
leads to the desired relations 



(Xyfy+X z f z )cly dz 


(5&) 







ds 


(5b) 


It remains to determine an analogous expression for the volume V 
between Z 0 and the body itself. The starting relation is Green’s 
theorem written in the form 


JJJ pCU^+uJdx dy dz 

T 

= -JJ P x ^ (U^x+cp)^ -JJJ x{ A [pCu^-HfJjj.) ] (pqy) + (pq> z )j- dx dy dz 


where r is the three-dimensional region enclosed by the surfaces £ 0 > 

and Z 2 . The factor within the braces, in the final term, is, however, 
zero throughout the region for steady-state, continuous, compressible flow 
conditions. To first order, the relation becomes 


JF 


P u dx dy dz 



§2 as 

OV 



2 

"P U + V 


cos (y>y) + w 

cos(v,x) 


cos(v,z) 
cos(v,x) J 


cos(v,x)d£ 


If, in the first term, the x-wise integration is performed and, in the 
second term, the boundary conditions on £ 0 are introduced and an x-wise 
integration by parts is carried out, one then gets 



dy dz = J fN(f,s)ds - U ro V + JJ f (Xyfy+X z f z )dy dz 


S 
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where N = N(x,s) is. a measure of the deviation of the body from the con- 
trol surface, that is, the distance measured normal to 2 0 between Zq 
and the body. Thus N = N(f,s) is the deviation of the body from the 
control surface at the base. From Green* s theorem, the last term is 


JJ fCVs+Vz)^ ^ = - Ij^x^dydz-i J X 


dff. 

8n 


ds 


and the volume formula becomes 


v - ^ff^l *• - I P £ (A, a )Jdy 

Co s C 2 


ds 


(6) 


Further simplification of equation (6) occurs if the trailing edge 
of the body lies in a plane perpendicular to the x axis. The first 
integral in the right member then becomes l times base area where l 
is body length and the term can then be re-expressed by means of equa- 
tion (5b). For trailing edges normal to the stream, volume then becomes 


V = 


ad/*' 


(i- f)‘ 


+ 1 r 




2U C 


;/ x s 

02 


ds 


(7) 


Solution of Variational Problem 


By means of the previous relations, wave drag can be minimi zed for 
prescribed values of incremental base area. A, and volume , V, external 
to Z Q . Since equation (7) is to be used, results involving volume V 
will be applicable only to configurations with trailing edges normal to 
the stream direction but this restriction need not be invoked when condi- 
tions are prescribed solely for the area A. The formal solution of this 
isoperimetric problem is achieved through the minimization of the 
expression 
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I = D - XA + UV 


( 8 ) 


•where D, A, and V are given, respectively, ty equations (tb) , (5b), 
ftnri (7) ari d are Lagrangian multipliers. Let X he replaced hy 

X + aT)(y,z) where r)(y, z) represents a variation that satisfies the same 
"boundary conditions as X on Gi and Cg and a is a constant. The 
necessary co ndi tion for minimization is that the expression 


jM* + 


X 

PcA 


f + 


2P, 


. U 

coco 


(l-f) 


, 1 Q 2 

+ 2 P 


(y‘ 


2 2\ 
' +Z ) 


} 


<3y dz + 



X 

P<xWx> 


f + 


M; 

2 PooUoo 




ds 


= 0 


holds for all 7. The following differential equation and houndary con- 
ditions must therefore he satisfied 


v 2 jx + 


A f + ^ 


Pco U oo 2p o3 U oo L 


(i-tf + 1 p 2 (Az 2 ) 


]} 




= 0 in S 


s[ X+ o 


X 


. u 

00 CO 


2poo^o 0 


= 0 on C, 


> (9) 


X = 0 on Cj 


J 


Equations (9) are closely analogous to conditions arising in the 
minimization of wave plus vortex drag for wings having given lift and 
given center of pressure. It should he remarked that the particular 
form of the term involving (y 2 + z 2 ) in the bracketed term of the har- 
monic differential equation is somewhat arbitrary and that, for example, 
the equation could he changed to one of the Poisson type with the term 
in question appearing in the right member as a constant times the 
Lagrangian multiplier p. 


For purposes of direct solution it is sometimes preferable to 
introduce the function ft(y,z) in equation (9) where 


ft(y,z) = x + 


x 

PxWx> 


f + 


2P3qUdO 



(10) 
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The equations to "be solved are, then. 


V 2 ^ = 0 in S 

ft = r~~ TT (y^z 2 ) on C 2 


Sn ^pooUoo ^ 


> (H) 


a = — ^— f + 


PboUcc 


2pooU ( 


- — I "(z-f) 2 + ^ ^(yS+z 2 )] on Cl 

coU ra L 2 J 


The function X thus differs from the harmonic function ft hy an amount 
given explicitly in terms of the characteristic surface S 2 . 

An Important expression for the wave drag can he derived from equa- 
tions (9) without further knowledge of the solution. For, from equations 
(4b) and (9), one has 


D = t If f + 2^k [ {l ~ t)Z + i p 2 ^ +z2 )]}^ ^ + 




C 2 


L Poo U co 


f + 


2pfcoU a 


U~f) 2 


ds 


and, by comparison with equations (5b) and (7), 

D = | A - | V 


( 12 ) 


In order to eliminate the Lagrangian multipliers, the solutions of equa- 
tions ( 9 ) or (ll) are needed. It is of interest to remark, however, that 
comparisons with results based upon combined-flow considerations show 
that for given base area and zero volume A can be identified with 
pressure in the combined-flow field and that for given volume and zero 
base area (closure) p Is the pressure gradient in the combined-flow 
field (see ref. 2). 


APPLICATIONS 


Equations (9) or (ll) determ in e the perturbation potential on the 
rear Mach envelope and thus lead to the evaluation of wave drag in equa- 
tions (4) . Particular applications are studied in this section. First, 
bodies with base area prescribed along a trailing edge are considered 
alone. Second, more particular cases of given volume and zero base area 
are treated. For given base area, it is proper to think of the body as 
extending downstream from the trailing edge along a cylindrical surface 
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having elements parallel to the stream direction. A semi-infinite shape 
then results. If the trailing edge is of supersonic type, all influences 
downstream of this edge can have no effect on the forward portion of the 
hody and, except for the hase drag arising from the unknown pressure at 
the trailing edge, the wave-drag calculations apply without regard to the 
cylindrical afterbody. 

The two final "bodies given here show how, for particular examples, 
the actual surface shape can be deter min ed from a knowledge of the func- 
tion X. Bodies of revolution are chosen for the e xa mples since the 
integral equations that appear in the analysis involve a single integra- 
tion. In general, double-integral equations will result and determina- 
tion of surface shape becomes correspon ding ly more difficult, if possible 
at all. 


Wings and Bodies Having Given Base Area 


The constraint on volume is now removed and the dependence of drag 
on base area alone is considered. From equation (10) , with p = 0, and 
equation (4a) , drag is 

D = ~ jj (% 2 +W z 2 )<3y <lz - —■ JJ ( W f zflz )dy dz + 

S S 


' t&J If (f ^ 2)ay 32 (13) 

S 

Green’s theorem can again be used to rewrite the first two integrals and 
the last term is simplified by virtue of the relation fy- 2 + f z 2 = (3 2 
In this manner one gets 


D 



GV 2 ^ dy dz 


S 



/ 

C3+C2 


f M ds + 

8n 


Pra 

2 


/ “H* 

C1+C2 


A 2 

Ho 


P 2 S 




dy dz + 


where, in the last term, q^^ is dynamic pressure 
tions (ll) , this relation reduces to 




From equa- 
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D = -A- f M f ds + -2f_ 0 2 S (1^) 

2U ro J ^ H, 

'■'l 

Further reduction of equation (l4) depends on a knowledge of the 
explicit solution of equations (ll) . A large and particularly interest- 
ing class of wings and bodies for which the solution is immediate is 
characterized by the condition f = const, on Ci. Since x = f(y,z) 
is the equation for the rear Mach envelope, the imposed condition implies 
that the outer curve Ti of this envelope lies in a plane normal to the 
stream direction. For example, all wings 'With plan forms having fore and 
aft symmetry satisfy this requirement as do also al 1 pointed configura- 
tions with subsonic edges so long as the nose and tail vertices deter min e 
a line parallel to the stream direction. The solution of equations (ll) 
is then 


ft = const. 


and equation (it) reduces to 


D = Jf- 

Ho 


(15) 


By means of equations (12) and (15 ) , A can be eliminated and drag 
expressed in terms of the geometry of the body or wing . The result is 


D _ A 2 

<3oo (H» 2 -DS 


(16) 


The simplicity of equation (l6) is remarkable and examples of its 
rather diverse applicability are given below. Before proceeding to these 
applications, however, it should be noted that a similar result applies 
to all planar wings for which the surfaces S x and S 2 intersect along 
any plane parallel to the z axis. The former condition that the curve 
Ti lies in a plane normal to the stream is thus relaxed so that 
x = f = m(y + b Q ) on the curve Ci, where m(< 0 ) is the slope of the 
plane of Ti relative to the stream direction. The solution of equa- 
tions (ll) is 


ft = (y+b Q ) in S 

Pco u oo 

and it follows that 

cos(n,y) on Ci 

Sn PcxMdo 
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After substitution of these results into equation (l4) , one gets 


D = (p 2 - m 2 )S 

4 <Jco 


(17) 


In this latter case, therefore, equations (17) and (12) give 


D_ = -A 2 = A 2 

<3oo S(p 2 -m 2 ) (tto 2 -l-m 2 )S 


( 18 ) 



Sketch (c) 


Ducted body of revolution -with pre - 
scribed end diameters .- This problem was 
first considered by Parker in reference 5 . 

As shown in sketch ( c) , a shape with ml nlmum. 
external wave drag is constructed so as to 
have an initial radius R^ and a final 
radius R x . In order that the previous 
linear theory should apply, the restriction 
is made that the ratio PlRx-Rgl/Z should 
be a small quantity. If the origin of axes 
is In the front face of the body, the fore 
and rear Mach surfaces are 


x = P(r-Ba) x-Z = -p(r-Rj 


. and the curve Cj. is a circle of radius R 0 where 

Hq = (Z+f3R 1+ pR 2 )/2{3 

From equation (l6) drag is 

p_ _ MRi 2 - R 2 a ) a * 

^ (Z+PRx+PRa ) 2 -^^! 2 


( 19 a) 


Equation (19a) is of particular interest since it represents a whole 
spectrum of results that extends from slender-body theory, for 3Ri/Z 
and pR 2 /Z small, to two-dimensional theory, for PRx/Z and large. 

The slender-body result leads directly to the familiar Karman ogive 
formula (ref. 6 ) 
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Elliptic plan form -with after - 
tody .- The problem of given base 
area along the rear edge of an ellip- 
tic wing was considered first by 
R. T. Jones (ref. 3)- The figure is 
a semi-infinite body with a cylin- 
drical shape drawn downstream of 
the rear edge of an ellipse, see 
sketch (d) . The equation of the 
plan form outline is assumed to be 



Sketch (d) 


a a U 2 


1 


and the enveloping Mach surfaces are determined completely by the fore 
and aft Mach cones with vertices along the supersonic-edged portion of 
the plan form, that is, the abscissas of the vertices lie within the 
region 


jx| < a 2 / ( a 2 +b 2 p 2 ) 1 /2 
The curve Ci has the equation 

2 o 

Z + , ! 

[(a (a/ P ) 2 

and is an ellipse with foci at (±b,0) . Equation (l6) then yields 


^ jta(a 2 +b 2 (3 2 ) 1/2 


(20a) 


If drag coefficient Op is based on plan-form area, equation (20a) can 
be re-expressed as 


°D ~ 


f-L ) 2 gl_ . , 

[p 2 +(V^) 2 f /2 


(20b) 


since the aspect ratio of the wing is AR = (Vb)/(jta). Perhaps the most 
convenient formula for comparison follows from equations (l9b) and (20a) 
if the drag of the wing is expressed in terms of the drag of a Karman 
ogive with the same length and base area. The ratio is given by 


_D_ = a 

( a 2 +b 2 p2) l/ 2 


1 


X/2 


[l+UPAR/4) 2 ] 


( 21 ) 
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The wave drag of the elliptic wing with cylindrical afterbody, in the 
limi t as aspect ratio approaches zero, is equal to the drag of the Karmen 
ogive and afterbody. For finite values of aspect ratio, the wave drag 
of the flat wing is smaller than that of the body of revolution, the 
initial deviation of the ratio from unity being proportional to (p®) 2 . 



I* — 2a — H 


Trapezoidal plan form with afterbody . - 
As a third example, consider a trapezoidal 
plan form of arbitrary taper ratio with 
base area along its trailing edge fixed. 
When the trailing edge is of subsonic type 
a cylindrical afterbody is assumed added. 

As shown in sketch (e), root chord is 
equal to 2a and span equal to 2b. The 
tip chord is 2d so that taper ratio Aq 
and aspect ratio ® may be introduced 
in the form 


Sketch (e) 


A 0 = d/a 



Sketch (f) 


JR = 2b/[a(l+A 0 )] 

So long as the leading edge of the 
wing is supersonic the characteristic 
trace Cx is as shown in sketch (f) and 
is composed of arcs of circles and straight 
lines, the radii of the inner and outer 
circles being a/p and d/p, respectively. 
The distance between the centers of the 
two outer circles is 2b. Once the lead- 
ing edge of the plan form is subsonic, the 
central circle of sketch (f) blankets the 
other parts of the figure and Cx is the 
circle of radius a/p. 


“by 


The area S is the sum of elementary geometric areas and is given 



(1+Ap) 2 

2 


PAl'cos a + a(l-A 0 2 ) 


where a, shown in the sketch, is given by 
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a = 


arc sin 


it 

2 


2(1-Ap) 

PAR(1+A q ) 


2(l-\>) < PAR(l+7^) 

2(l-\>) £ (l+\)) 


From equations (l6) and (l9h ) , the min i mum drag of the trapezoidal wing 
relative to the drag of the Karman ogive of equal length and base area 
is 


_D^ 

Dk 


*+(l+\,) [(l+A 0 ) 2 3 2 iR 2 -4(l-A 0 ) 2 ] 1/2 -2(l-?\ 0 2 ) are cos \ 2( ^P } ■ 


when 

2(l-7so) < pA* (l+*o) 



when 

2(i"\>) ^ m (1+Ao) 


Special cases of interest are: 
Rectangular wing (A,-, = 1) 



= 1 
Ife 1 + (4pA* )/jt 


(23) 


Diamond wing (A q = 0) 


J3_ 

% 


1 

1 + i [( P £*) 2 - 4] 1/2 



2 £ piR 

> 


(24) 


= 1 , PiR ^ 2 I 
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Sketch (g) shows a plot of D/Dg against PAR. for the elliptic, 
rectangular, trapezoidal, and diamond plan forms. Base area and length 
of the wings are equal to these parameters for the Karman ogive. For 
large values of PAR the relative drag decreases as 1/ pAR . As the wing 
plan forms become slender, drag of the elliptic and rectangular wings 
approaches in the limit of vanishing PAR the drag of the ogive. The 
tapered wing, on the other hand, has a value of drag equal to that of the 
ogive for all values of taper ratio and aspect ratio satisfying the 
ineq uali ty 2(l->o)= PAR (l+Ao) * This relation is satisfied so long as 
the edges of the wing are subsonic. Changes in sweep angle of the lead- 
ing and trailing edges produce no further change in the minimum drag of 
the configuration so long as the base area is held fixed. The value Dr 
is the ml nimum drag for all such configurations lying within the fore and 
aft Mach .cones from the nose and tail of the w ing . Except in the case of 
the rectangular plan form, the curves of D/Dg have zero slope at their 
peak values. 



Sketch (g) 
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It should be noted that among the optimum configurations considered 
here the Karman ogive has the greatest value of wave drag for given 
length, base area, and Mach number. The generality of this observation 
is apparent from equation (l 6) . Let the family of wings include all plan 
forms for which the curve Tx lies in a plane normal to the stream direc- 
tion; the maximum value of minimum drag occurs when S is a minimum and 
this occurs when the trace C 2 is the circular curve associated with the 
Karman ogive of equal length and base area. 


Yawed elliptic plan form with 
afterbody .- For given base area, the 
drag of an elliptic plan form at 
angle of yaw ij r can be calculated 
from equation (18) . In order to 
justify this statement it is suffi- 
cient to show that the characteristic 
curve Tx lies in the plane 
x = f = m(y + b Q ) . The trace of Tx 
in a yz plane is, in fact, another 
ellipse and the dimensional relation- 
ships between the plan form and the 
trace are as shorn in sketch (h) . It 
is convenient in the derivation of 
these results to proceed inversely 
and to determine the plan form as an 
envelope of curves given by the inter- 
sections in the xy plane of fore and 
rear facing cones with vertices on Tx. 
Since the streamwise position of the 
origin is of no direct significance, 
the space curve Tx may be assumed 
given by 



B > C 


x = my f m < 0 



1 



The Mach cones with vertices at the 
point (xx, yx, z x ) on Lx are 


(x-xx) 2 = p 2 [(y-yx ) 2 + (z-zx) 2 ] ( 26 ) 


Sketch (h) 


where 

zx 2 = (B 2 -y x 2 ) 

B 2 


xx = myi , 
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The parametric equation of the envelope in the yz plane is given by 
equation ( 26 ) with z = 0 and the y x derivative of the same expression, 
that is, the relations 

B 2 [(x-my 1 ) 2 - P 2 (y-yi) 2 ] = p 2 C 2 (B 2 - y:i 2 ) 

[B 2 m 2 -p 2 (B 2 -C 2 ) ]yi = (mx-p 2 y)B 2 
The envelope is, therefore, 

(B^C^x 2 2mB 2 xy + (m^-^C^y 2 = C 2 [3 2 (B 2 -C3 2 ) - m^ 2 ] (27) 

Equation (27) represents an ellipse so long as the initially chosen m 
satisfies the inequality 

m 2 < p 2 (B 2 -C 2 )/B 2 (28) 

The elliptical plan foim is fixed by its major and minor axes and 
angle of yaw. The relationship between the plan form and trace curve is 
more conveniently carried out, however, in terms of the three quantities 
Z, b, p where Z is streamwise length of the plan form, 2b its width, 
and x = py is the line passing through the points on the plan form 
where y = ±b. 

Elementary calculations performed with equation ( 27 ) yield the 
following relations 


Z/2 = (BSmS+pec 2 ) 372 

( 29 a) 

b = (B 2 ^ 2 ) 372 

( 29 b) 

mB 2 

11 ~ B 2 — C 2 

( 29 c) 

-2mB 2 

tan 2t}t = — — 

(B 2 -C 2 ) - (m^ 2 -^ 2 ) 

( 29 d) 


In sketch (h) the plan form is also circumscribed by a parallelogram 
with sides inclined at the Mach angle. The equations of these lines are 

x = Py ± (g-m)B x = -Py ± (p+m)B 

from which it follows that their outermost intersection points are at 
y = ±B and the line connecting the intersection points is x = my. 
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The above results thus show that the Mach lines circumscribing the 
plan form can be used to determine the span of the trace of Tx and 
the angle of inclination of the plane of fx. The span of the plan form 
is, moreover, equal to the distance between the foci of the elliptic 
trace . 

From equations (18) and (19b) the drag of the wing and afterbody 
relative to the drag of the Karman ogive of equal length and base area 
is given by 


p _ I s - _ B 2 m g +p 2 C 2 
4BC(p 2 -m 2 ) BC(0 2 -m2) 


The results can be summarized as follows 


_d_ = i4B 2 i 2 +[(i+p 2 ^ 2 ) 2 "^ 2 p 2 |^3 1/2 1/2 1 

l“P 2 i 2 +[ (1 +p 2 £ 2 ) 2 -4p 2 !i 2 £ 4 ] 1/2 [ (l+p 2 | 2 ) 2 -4p 2 p 2 £ 4 ] 1/2 

(31) 


m = 


*c(ln a S a ) Va ' 


I = 2b /l 


tan 2\|r 


2p| 2 

l-E 2 


(33) 


Sketch (i) is a plot of D/D g 
against angle of yaw for M = */2 
and Al = 4/*, 4, 8. The smallest 
of these values of aspect ratio 
corresponds to a circular plan 
form and obviously must be inde- 
pendent of i|r; the drag ratio is 
D/Dk = */2/2 and this is in agree- 
ment with equation (21) for the 
special case of the circular wing. 
Several limiting forms of equa- 
tion ( 31 ) are of interest in show- 
ing the variation of drag. For 
example, when p = 0, the plan 
form is unyawed, JR =4 £/it and 



Sketch (i) 
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D _ 1 1 

(l+£ 2 g 2 ) 1/2 [l+(:rtpiR/4) 2 ] 1/2 

as given in equation (21) . This relation furnishes the values of D/Dj£ 
in sketch (i) at i|r = 90°. If pg = 1 and £ ^ 0 one has M = co 
ijr = arc tan g and 

_D_ _ 1 

^ (l-p 2 tan 2 i|r) 3/ 4 

This is the general drag relation for the yawed wing when aspect ratio 
becomes infinite. It is to he noted that drag r emain s f ini te except when 
the angle of yaw is equal to the free- stream Mach angle. In sketch (i) 
the drag curve for infinite aspect ratio must therefore have a singularity 
at ljr = 45°. If pg 1 and g 0 so that aspect ratio remains finite, 
it follows that i]r 0 and D/Dk = 1. 


Ducted Body of Revolution of Given Volume 




Sketch (j) 


External wave drag of a quasi- cylindrical 
body of revolution is to be minimized. Bee 
sketch (3) . The radii of the body at the 
nose and tail are assumed equal (=R) and the 
volume, V, of the body external to the refer- 
ence cylinder is a known constant. The value 
of the minimum drag and, as will be developed 
in the next section, the evaluation of the 
body shape can be determined from the solu- 
tion of equations (ll) . The equations of the 
fore and rear Mach envelopes Eire given 
respectively by 


x = p(r-R) x = f = l - P(r-R) = p(2h+R-r) (34) 


where the constant h (=z/2f3) is introduced for convenience. The curve Pi 
is a circle of radius r = R + h in the plane x = {3h and equations (ll) 
thus take the form 



HACA EH 3289 


21 


V 2 *! = 0 


1(1=^ 


Sr 2 Pco U a 


AP 

a = — — h + 


for r = R + h 
at r = R 


( 35 ) 


nP £ 


Pco U a 


SPcoUa 


[ 


h 2 ^ (*+*)' 

h + 


at r = R + h 


The solution of equations (35) can he expressed in the form 


ft = hi + k2ln[r/(R+h) ] 
where, from the boundary conditions. 


ki = 


AP pp 

h + 


PcoU c 


2pfcoUco L 


2 (R+h) ‘ 

h + 


k2 = 


|ip 2 R g 

2 PcoUoo 


This result, together with equation (10) , determines the function X 
once the relationship between A and p is established from the body 
closure condition. If the known quantities X and f are substituted 
into equation (5b) and A set equal to zero, the desired relation is 
seen to be 

A = hpP (36) 


and the final form for ft is 


ft = 


HP* 


2 PtoU a 


3h 2 +i?l h . ) : 


• + R 2 !n 


R+h 


( 37 ) 


Drag can now be evaluated directly. A parallelism with the deriva- 
tion following equation (13) can be maintained if equation (4a) is used 
and the potential is written 


X = fl(y,z) - F(y,z) 


( 38 ) 


where, from equation (10), 


= — f + — H: — 

Poo u co 2 Poo U co 



+ \ P 2 (y 2 +z 2 ) 


( 39 ) 
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Drag is given by 


D = ' 2 " JJ (% 2+fi z 2 )^y ~ P c a JJ (Fyfly+F z n z )<3y dz +^° JJ ( F y 2 +F z 2 )dy dz 

S S S 


If Green’s theorem is used in the first two integral terms , along with the 
boundary conditions , one has 


D = ~ JJ (F y 2 +F z 2 )dy dz 

S 



+ 



pP 2 _S_ (y 2 +z 2 ) 

2pJJ m ^ 2 




Equation (1*0) applies in general. For the particular problem of the 
ducted cylinder, equations ( 3 !*-) and ( 39 ) give 


F = 


nP 2 

2Poo U co 


3 h 2 + (R+h-r ) 2 



and substitution in equation (to) yields 


D = | tUL. | ^ (2Rfh)[(R+h ) 2 - 3R 2 ] + R 4 Zn ^ j (1*1) 

The multiplier p can now be eliminated between equation (12), with 
A = 0, and equation (1*1) . The final result assumes a relatively concise 
form when written as follows 


_D_ = V 2 
%> Z 4 C(n) 


(tea) 
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where 


C(cr) = ~ ^(l+lkr)[(l+2ff) 2 - 12<j 2 ] + 6kaS n 

As in the case of the quasi- cylinder with given base area, these 
results cover the entire spectrum of fineness ratios and yield, in their 
limi ting forms, the results of two-d im ensional airfoil theory (biconvex 
section) and slender-body theory. The latter case, which is the Sears- 
Haack slender body (refs. 7 and 8), corresponds to a = 0. Equation (42a) 
then becomes 

P S-H _ 128V 2 
%> Til 4 

The above problem was considered previously by Heaslet and Puller 
(ref. 9) without recourse to the present techniques but, rather, after 
expressing drag in terms of the source distribution that could be assumed 
to generate the external shape of the body. In this approach it becomes 
necessary to find first the source -distribution function, under minimizing 
conditions, and to calculate drag and volume subsequently. The details of 
the calculations are thus less direct since the desired quantities are 
expressed as integrals involving the hyperbolic influence function of the 
supersonic source. In reference 9, the function C(cr) of equation (42a) 
appeared in the form (in a slightly modified notation) 

C(cr) = |J" 1 [(TH-2a)(l-ri + 2c)] l/2 [T 1 (l-Ti)E-a(l-4a)(K-E)]dT 1 (43) 

o 

where K and E are elliptic integrals of the first and second kind, 
respectively, of modulus 


(42b) 



k 


-(t)+2u) (l-r)+2cr) 


, 1/2 


The immediate advantage of equation (42a) is, of course, the natural 
advantage provided by any analytic representation with its precise deter- 
mination of magnitude and rate of change. From a disparate point of view, 
the equivalence of the two results gives not only a new fundamental iden- 
tity in the theory of elliptic functions but also indicates a method 
whereby further identities can be generated. From the standpoint of 
direct application, however, the results of reference 9 remain unmodified. 
The calculations that were used to plot the variation of C(cr) were found 
to check to at least four significant figures with the present formula 
and thus provided a satisfying confirmation of the numerical techniques 
vised in the original evaluation. 
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Determination of Surface Shape 


In the calculation of minimum drag for the two families of ducted 
bodies, the value of the perturbation .potential on the rear Mach envelope 
becomes available in explicit form. It is therefore possible to seek the 
source distribution function B(x) generating the disturbance field on 
the rear surface; once B(x) is known the body shape can be calculated 
as the final step in the analysis of the mi nim um drag bodies. When the 
source distribution is confined to a single line segment, B(x) appears 
as the solution of a single-integral equation. The origin is conven- 
iently moved to the starting point of the source distribution (x = -PR 
of sketch (h)), so that B(0) -= 0. For bodies of revolution, the 
perturbation potential is known to be given by 



B(x 1 )dx 1 


(VO 


Since B(0) = 0, differentiation and integration by parts yields 


x-Pr r -, 1/2 

/ [^|] 


-S-Px + Hr " 5^ 


where the prime denotes differentiation. The rear surface is 


x = f = l + 2PR - f3r 


and on this surface 


*r = -0<Px + q>r 


Hence, the integral equation 


^ SR - aPr " l+2(3R-X! 


v i — - r -*■ r - 

2jrrX r = J 


I+2PR-2Pr-x 1 J 


B' (x 1 )dx 1 


(^ 5 ) 


results where 2irrX r can be assumed known. The transformations 
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t = 1+ 20R-2gr 


G(x) = (z+2PR-x) l/a B' (x) 


H(t) = 2iccX r (r) 


m 


lead to the Abel type equation 


=(t> -/ 


G(x 1 )dx 1 


J o (t-x x ) 


1/2 


( 47 ) 


and its solution is known to be 


G(t)=ilf 

' X dtj 


H(t 1 )dt 1 
1/2 


o (t-t x ) 


( 48 ) 


For the case of given base area (p = 0), the solution of equations (ll) 
was ft = const., hence, from equation (10) 


X = - 


PooU* 


f + const. 


and 


X r = 


_AP_ 

P«Hx 


From equations (12) and (15) it follows that 


P 2 S 


so that one has 

H(t) = (l+20R-t) 

p 2 S 
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Substitution into equation (48) yields G(t) and from equations (46) 
the source distribution functions (with x as in sketch (h) ) 


b'M-M 


7 - 2x 


P s [ (x+§R) (7+gR-x ) ] 1/2 

B(x) = ^ [(x + gR)(7+gR-x )] l/2 
g^S 


m 


For the case of given volume, equations (37 ) > (38), and (39) yield 


rX-p = 




2PooU 0 


[R 2 + 2(R+h)r - 3T 8 ] 


and the determination of the source distribution function follows directly. 
The final form of this function is 


B(x) = U “ V 


27 4 C(a) 


-f( 7 - 2 x)[ (x+gR) (l+gR-x) ] 1/2 - 2 g 2 R 2 arc cos X ( 50 ) 

t 7+2 gR J 


The geometries of the two families of bodies were calculated and 
shown in reference 9 and trill not be repeated here. Comparable results 
to those of sketch (g) are, however, given in sketch (k) . For the case 
of given volume, the drag of the quasi -cylinders relative to that of the 
Sears-Haack slender body is, from equations (42a) and (42b), 



Sketch (k) 


D it 

D s _ h ” 128C (cr) 


(51) 


For the case of given base area the 
drag formula is reduced to the form 
it would take in reference surface 
theory. Under the assumption that 
R^ deviates slightly from R x , 

equations ( 19 a) and ( 19 b) yield 


_D_ _ 1 

Dg; 1 + 4d 
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CONCLUDING REMARKS 


As given "by equation (4a) , drag is equal to the kinetic energy in a 
two-dimensional flow field of unit thickness. Intuitively, the results 
are r emin dful of the similar expression for the vortex drag of a lifting 
wing in subsonic flow theory. In this latter Trefftz plane analysis, the 
optimum field is represented by a harmonic function in the entire cross- 
flow plane external to the trace of the body whereas, in supersonic flow, 
an additional, .outer, bounding curve Cj. appears and affects the induced 
velocities and the wave drag. The flow field associated with minimum drag 
due to lift remains a harmonic field on the rear Mach envelope in super- 
sonic flow (ref. 1) and, as shown in the present paper, the flow in the 
nonlifting case is also a harmonic field with known additional effects. 
Another point of similarity between the subsonic and supersonic theories 
also appears: the boundary conditions determining the induced field are 

given along curves that are uniquely determined by the wing or body but 
the details of plan-form shape are not determined uniquely from the 
bounding conditions themselves. Thus, in Prandtl’s and Munk's vortex 
theory, the flow in the Trefftz plane can be calculated but the chordid.se 
dimension of the wing disappears entirely; all wings with equal span have 
the same ml nl mum drag and the same spanwise distribution of loading for 
minimum drag. It is obviously impossible to solve for details of the 
chordwise load distribution in this case. In supersonic theory the condi- 
tions for minimum drag are given on a surface displaced from the body and 
a loss in knowledge as to plan-form shape is thus incurred. In the case 
of the trapezoidal plan forms with given base area, minimum drag and 
induced velocities on the Mach envelope remain fixed for the entire f ami ly 
of wings with subsonic leading and trailing edges. The solution of the 
problem allowed for no storage of information about taper ratio and aspect 
ratio once the wing was swept behind the Mach cone from the vertex. There 
remains the final question, however, as to whether a multiplicity of solu- 
tions can be found. Suppose a planar ling exists and consider the shape 
equal to the difference between the wing and the Karman ogive. From equa- 
tion (l4) the drag of the latter shape is zero. But for nonvanishing 
thickness, wave drag can never be zero. The logical contradiction implies 
that the Karman ogive is the only minimum drag configuration within the 
family of trapezoidal plan forms with subsonic edges. 

Following equation (12) it was remarked that the parameter A could 
be identified "with pressure in the combined flow field. An interesting 
interpretation of the value of this parameter follows from the use of 
equation ( 15 ) . Eliminating drag between the two relations, one gets 


2qco " p 2 s 


( 53 ) 
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The left-hand expression is half of pressure coefficient in the combined 
flow field. Equation ( 53 ) states that the value of -~h / necessary to 

calculate minimum drag is equal to the pressure coefficient predicted by 
linearized unidimensional flow theory in a duct bounded internally and 
externally by the characteristic traces, that is, the curves C x and C 2 
of sketch (b) . 


Ames Aeronautical Laboratory 
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Moffett Field, Calif., June 2k, 1957 
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